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Experimental nuclear moments of inertia at high spins along the yrast line have been determined 
systematically and found to differ from the rigid-body values. The difference is attributed to shell 
effects and these have been calculated microscopically. The data and quantal calculations are in- 
terpreted by means of the semiclassical Periodic Orbit Theory. From this new perspective, features 
in the moments of inertia as a function of neutron number and spin, as well as their relation to 
the shell energies can be understood. Gross shell effects persist up to the highest angular momenta 
' observed. 
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> i I. INTRODUCTION 

o . 

, Microscopic calculations based on the rotating mean field (self-consistent cranking model) describe the rotational 
' energies well. Although there are numerous detailed comparisons between experiment and calculations in the liter- 
, ature, a qualitative understanding of the rotational response has not yet been reached. Two idealized models have 
been studied in detail: the Fermi gas and the harmonic oscillator P]. The Fermi gas does not take into account the 
" ' ] shell structure. Its yrast line is characterized by Jrig, the moment of inertia of a rigidly rotating mass distribution 
' (the rigid-body moment of inertia). The harmonic oscillator has a very specific shell structure, which is different 
from the one of the nucleus. At moderate angular momentum and equilibrium deformation it has also the rigid-body 
moment of inertia. However, we know from experiment that the moment of inertia along the yrast line is in general 
T-H different from jTrig, indicating that there are currents in the rotating frame. The moments of inertia observed at low 
. spin reach only about 40% of the rigid-body value. This reduction has been commonly attributed to the presence of 
• pair correlations Q . If the nucleus were a macroscopic superfluid it would have the irrotational flow pattern, which 
corresponds to strong counter currents in the rotational frame and moments of inertia J^irrot that are substantially 
smaller than J^rig- The observed moments of inertia lie in between J^irrot and J^rig- This can be explained by the 
fact that the coherence length of nuclear pairing is larger than the nuclear size, which prevents the nucleus from fully 
developing the irrotational flow pattern of a macroscopic superfluid. 
^ ' However, the moments of inertia observed at high spin also substantially deviate from the rigid-body value, although 
(~| the pair correlations are expected to be small. This observation is corroborated by self-consistent cranking calculations 
that assume zero pairing (see e. g. which find strong deviations from Jrig in accordance with experiment. 

Nevertheless, the question of long standing interest in nuclear physics, how is angular momentum generated in nuclei, 
is not yet answered in a satisfactory way. Shell effects (i.e., the bunching of single-particle energy levels) play an 
?H important role at high angular momentum. Superdeformation ^5j, band termination Q and uniform rotation about 
a tilted axis are dramatic manifestations of shell effects. The shell effects determine the moments of inertia along 
the yrast line, the importance of which was pointed out by Pashkevich and Frauendorf long ago. The deviations 
from the rigid-body moment of inertia imply that the flow pattern must substantially deviate from the current of a 
rigidly rotating mass distribution, i. e. there are strong net currents in the frame rotating with the density of the 
nucleus. In the present paper we investigate the deviations of the moments of inertia from Jrig^n. a systematic way 
and give a qualitative interpretation in the frame of the semiclassical Periodic Orbit Theory which relates the 
shell structure to classical periodic orbits in the nuclear potential. This approach becomes a very instructive tool for 
interpretation if one is interested in the gross structure, which is determined by only the few shortest orbits. We take 
a fresh look at the angular momentum generation in nuclei through this perspective. 

In section^ we extract the experimental moments of inertia at high angular momentum, where the pair correlations 
are expected to be negligible or at least drastically reduced and calculate their deviations from the rigid-body value. 
We find substantial deviations from the rigid-body moments of inertia. In section ITTll these experimental deviations 
are compared with mean-field calculations, which assume zero pairing and employ the shell-correction version of 
the self-consistent cranking model (Strutinsky-type calculations in a rotating Woods-Saxon potential). We find that 
the calculated deviations from the rigid-body moments of inertia are rather similar to the experimental ones. We 
also extract experimental shell energies which are compared with the calculated ones. In section Hvl we review the 
elements of the Periodic Orbit Theory which we need to explain some features of the gross shell structure at high 
angular momentum . We find that the shell energies and the shell moments of inertia are correlated and determined 
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by the contributions from the shortest classical periodic orbits. The competition and interference between the orbits 
in the meridian planes and the equatorial plane of the deformed potential account for the most conspicuous features 
of the shell energies and of the deviations of the moments of inertia from the rigid-body value. They also explain 
the appearance of regions of high-spin isomers. The deviations of the moments of inertia from the rigid-body value 
indicate the presence of currents in addition to the simple flow pattern of a rigidly rotating mass distribution, i.e., of 
net currents in the rotating body- fixed frame of reference. A first report of our results was given in ref. 0. 



II. EXTRACTION OF MOMENTS OF INERTIA AND SHELL ENERGIES FROM THE DATA 

We are interested in the gross shell structure, not in the details of the dependence of the energy on the particle 
numbers Z and A'^ and the angular momentum /. Accordingly, we characterize the yrast line (sequence of states 
with minimal energy at a given I) by only one parameter, the moment of inertia J', assuming for the energy the 
/-dependence. 

The moment of inertia we want to determine for each nucleus is the yrast moment of inertia in the terminology 
introduced by Bohr and Mottelson (see Fig.l in It characterizes the increase in energy of the yrast levels with 

spin / and represents the global characteristics of the generation of angular momentum in the nucleus. It should be 
emphasized that the yrast line does not necessarily coincide with one rotational band. It may be composed of pieces 
of crossing bands or not contain bands at all. Since the expression ^ describes well the average / dependence of the 
yrast at high spin (see below) , there is no need to distinguish between kinematical and dynamical moments of inertia. 
In Fig. n] examples of the yrast energies as functions of /(/ + 1) are given for a "rotational" nucleus (top, ^^'^Er) 
and for a "non-rotational" nucleus (bottom, ^^'-'Dy). While the plot is smooth for the rotational nucleus, due to 
the presence of long rotational bands, there are irregularities for the non-rotational nucleus, which reflect alignments 
of individual nucleons with different angular momentum. Besides these fluctuations, there is an approximate linear 
relationship at high spin in both cases. Since we are interested in the average behavior along the yrast line, we fit the 
linear relation to the yrast line at the highest spins, which will give the yrast moment of inertia. We do not fit 
separate bands, which would give a band moment of inertia. This procedure is applied to the rotational nuclei and to 
the non-rotational nuclei for which the yrast line is irregular. 



A. Quenching of the pair correlations at high spin 



Generating angular momentum destroys the pair correlations. This is in analogy to the transition from the paired 
to the normal phase of a superconductor in a magnetic field. In a finite nucleus, the transition is less distinct and 
stretched over a substantial interval of angular momentum. Hence, there is a question as to what extent one can extract 
from the high-spin data moments of inertia that characterize the rotation of the unpaired nucleus. The quenc hing 
of pairing by rotation has been studied in a number of pubhcations. See e.g. the recent studies 0, ^2 0|j 
which refer to the extensive earlier work. The following picture emerges. The typical pairing effects, which include 
the large reduction of the moment of inertia at low spin, are caused by strong correlations in the occupation of a few 
levels near the Fermi surface. They are well accounted for by the mean-field (Hartree-Fock-Bogoljubov) description 
of pairing. These so-called static pair correlations are destroyed by breaking few pairs. The destruction of the static 
correlations is seen as a change of the slope of the yrast sequence, which typically appears in the range 14 < / < 26 
for the considered nuclei. On top of the static pairing there are the dynamic pair correlations, which involve the 
levels far from the Fermi surface. They decrease slowly with increasing angular momentum, remaining substantial 
in the experimentally accessible spin range. They reduce the moments of inertia by about 5%, where the reduction 
is not sensitive to the detailed structure of the states near the Fermi surface. This explains why unpaired cranking 
calculations (see e.g. Q) are so successful in describing the level structure in the yrast region at high spin. Hence, 
we try to extract the moments of inertia from the region where the static pairing is destroyed. These values are 
expected to be close to the moments of inertia of the unpaired nucleus, with a slight overall reduction. In particular, 
the experimental moments of inertia will show the shell structure of the unpaired nucleus. 
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B. The extraction method 



As seen in Fig. ^ the experimental excitation energies E{I) when plotted vs /(/ + 1) are roughly linear at high 
spin, where the static pairing is destroyed. The experimental curve bends sharply down at low spins (J < 15 and 
most often / < 10, especially in lighter nuclei). This is caused by the onset of the static pair correlations, which 
reduce the moment of inertia. The examples are typical for all the studied nuclei. Since we want to eliminate the 
pairing effects, we fit the yrast lines with a straight line for the highest lOh where possible, including both odd and 
even spins, as shown in Fig. ^ If negative parity levels enter the yrast line they are accepted. To avoid interference 
from low-spin pairing effects, we include only nuclei that are known up to sufficiently high spins, using the following 
criteria: the highest known spins should be greater than Ibh for nuclei with neutron number < 50, greater than 
16h for 50 < iV < 82, and greater than 17h for N > 82. The minimum spins entering the fit should be greater than 7h 
for N < 50, greater than lOh for 50 < iV < 82, and greater than 12h for N > 82. Generally the maximum spins are 
considerably higher than the limits above. The experimental yrast energies are taken from the ENSDF files using the 
Isotope Explorer program jl5| . for nuclei with even Z in the mass range from 80 to 220. From the straight line of the 
fit, E{I)fit — Eg + 1(1 + l)/2j7e2;p, we derive the experimental moment of inertia J^exp and the ground-state energy 
Eo of the unpaired nucleus. Fig. |21 shows the difference E{I) — Efit{I) for a representative selection of nuclei. At 
high spin, the curves largely flatten out near the zero line. There is a distinct drop of the yrast levels with decreasing 
/, which is interpreted as the gain in energy caused by the static pair correlations. The onset of this correlations is 
seen both in rotational and spherical nuclei. The difference Eo — E(0) may be interpreted as the experimental ground 
state pair correlation energy. It turns out to be rather different for the examples. It would be interesting to see if 
theory can reproduce these variations in a systematic way. The deviation of the upper portion of E{I) — Efit{I) from 
the zero line reflects the uncertainty of the experimental moments of inertia. 

The rigid-body moment of inertia is calculated using the formula 



where Pg is the ground-state deformation as calculated by MoUer et al. fl6j . Since the rigid-body moments of inertia 
are rather insensitive to the deformation, taking the ground-state deformation instead of the actual deformation at 
high spin will have only a small effect. We then calculate J^yrast — Jrigi the deviation from the rigid-body moment 
of inertia, and plot it as a function of neutron number. The experimental moments of inertia, and their deviation 
from rigid-body values for the nuclei with small or normal deformation are shown in Figs. |31andJ^respectively. The 
superdeformed bands of the A — 150 region are shown in Fig. O For J, we use the estimates of pTj . 
The shell energies at spin / are calculated in the Strutinsky sense [Tg by means of the formula 



where i? is the experimental binding energy and i?LD is the liquid-drop energy at the ground-state deformation taken 
from 16]. Fig. shows the ground-state shell energies and Fig. [7|the shell energies at / = 20. The energy Es/j also 
contains the fluctuations of the pair correlation energy since the smooth pair energy is included in ^ld- 

As seen in Fig. 01 the deviations of the moments of inertia from the rigid-body value are dramatic for small and 
normal deformation. Rough numbers for the experimental deviations are -25% for N ~ 75, -1-20% for A^ ^ 86, -15% 
for A^ ^ 95, -40% for A^ ^ 100, and -65% for A^ ~ 120. Note the positive value, which cannot possibly be attributed 
to pair correlations. We are going to interpret the deviations in terms of the shell structure. Here, we want to argue 
that they are not caused by the limited angular momentum reached in experiment. Fig. |S1 shows the maximal value 
of / of each nucleus used in the analysis. There is no obvious correlation with Fig. 01 except near the closed shells 
N = 82 and 126, where both the moments of inertia and the maximal values if / are small. In Fig. |51 we show only the 
nuclides for which the maximum spin is larger than 20. They represent approximately 65% of all nuclei considered. 
The overall picture of deviations of the moments of inertia from the rigid-body value is very much the same as in 
Fig. 01 Since for the subset of nuclei with the highest spins the pairing correlations are, in any case, expected to be 
weaker than at lower spin, the similarity of the structures in Figs. 01 and indicates that neither the variations in 
experimental spins nor pairing correlations are creating that structure. Hence we conclude that Figs. 01 and El show 
the contribution of the shell structure to the moment of inertia and are not related to pairing effects. 



J„g[h'^/MeV] = 0.01253^5/^(1 + + 0.048225A, 



(2) 




(3) 



III. COMPARISON WITH QUANTUM MECHANICAL CALCULATIONS 



In this section we demonstrate that the experimental shell-energies and moments of inertia can be accounted for 
by a quantal treatment of the rotating mean fleld without pair correlations. 



4 



A. Cranked Woods-Saxon-Strutinky calculations 

For a realistic quantum mechanical calculation we use the Strutinsky-type shell correction for rotating nuclei |l9l |. 
The single-particle levels are calculated at the rotational frequency uj by diagonalizing the cranked deformed Woods- 
Saxon potential 

h' = t + V{a, ai)-uj- j, (4) 

where j is the single- particle angular momentum . For the single-particle potential a simple generalization of the 
spherical Woods-Saxon potential is used [20I l2l| . It is supposed that the potential depends on the distance from the 
surface ^d(r). 

V{r)^Vo/[l + expildir)/a)], (5) 



Vo = -V,nean{l ± c'^^°^^^), (6) 

where the upper (lower) sign refers to protons (neutrons). The distance ld{r) is determined numerically. The sign of 
ld{r) is taken to be negative inside the nucleus and positive outside. In case of the sphere Idi^) = fspft — Ra, where 
i?o roAa. The spin-orbit interaction is expressed in terms of the central potential V as follows 

V^-^-^VV-lap]. (7) 

The parameters of the potential are given in tabled They are taken the same for all nuclei considered as the averaged 
of the four sets of the parameters given for the rare earth region in ref. . The parameter rp enters in the definition 
of the potential through the volume conservation condition. 

The axially symmetric shape is characterized in the meridian plane in the coordinate system in which one of the 
two families of the coordinate lines are the con- focal Cassini ovals . For small deformations the shapes are almost 
spheroidal and the parameters a and are close to the usual quadrupole and hexadecapole deformation parameters 
e and £4. The shell energy (Strutinsky shell correction) Esh[oi, 04, w) is obtained from the single-particle levels by the 
Strutinsky's averaging procedure [l^ . The deformation parameters are determined by minimizing with respect to a 
and a4 the total Routhian (energy in the rotating frame) 

E'{a,ai,uj) = Esh{.a,ai,u) + ELD{a,ai) - —Jrig,,,ia, a^), (8) 

at the given rotational frequency uj. We consider the two possibilities that the rotational axis uj is perpendicular to 
the symmetry axis {v = _L) and that it is parallel {v = ||). The angular momentum is the expectation value of its 
projection on the rotational axis, 

J^{uj) =< uj\j^\uj >, (9) 

where = _L or || and \uj > is the lowest configuration in the cranked Woods-Saxon potential The expression for 
the liquid-drop energy E^jj is given in p3 |. The calculated moment of inertia at the frequency ui is 

^,(^) ^ M^. (10) 

UJ 

In order to emphasize the role of shell effects, we call the deviation of the moment of inertia from the rigid-body value 
in absence of pairing, 

Jsh,iy{^^) = Ji/(w,a, 04) - Jrig,iyia,a4), (11) 

the shell moment of inertia. We calculate the shell moments of inertia for uj = 0.3 MeV/h in the range of nuclei for 
which there is experimental data. The nucleus is allowed to rotate around the symmetry axis and around an axis 
perpendicular to the symmetry axis, and the shape parameters are optimized separately for prolate and oblate shape. 
Figs. ^1 shows the four types of shell moments of inertia . Out of the four calculations, the mode with the lowest 
total energy is chosen to calculate the optimal shell energies and the optimal shell moments of inertia. The results 
are shown in Figs. II II and [T^ 
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B. Discussion 

The experimental shell moments of inertia are shown in Fig. 0] as functions of neutron number TV. The proton 
numbers Z are represented by the different symbols. The calculated shell moments of inertia are shown in Fig. 1111 In 
both cases the yrast moments of inertia deviate substantially from the rigid-body values. The variation of the shell 
moments of inertia with neutron number is very similar in the experiment and in the calculation: there are dips at 
the spherical neutron magic numbers 50, 82 and 126 and peaks just above and below N=82 then lower values around 
N=90 where the deformation sets in, then another peak around N=110 which is the region of high-K isomers. We will 
see later that these features can be related to particular properties of the nucleus (such as deformation, closed shells 
and axis of rotation). Overall, the experimental shell moments of inertia are shifted by about 10% to the negative side 
compared to the calculated ones. We attribute this shift to the residual pair correlations, which are mostly dynamical. 

The similarity between the moments of inertia extracted from the experiment and from the calculations that do 
not include pairing effects is a confirmation that the differences, J — Srig, are not due to pairing effects, but are the 
manifestations of shell effects. In fact, the present comparison is just another example of the observation that the 
details of the rotational response at high spin are well reproduced by calculations without pair correlations, which 
often give moments of inertia that are very different from the rigid-body value. One case is the moments of inertia 
in high-K multi quasiparticle bands Q. The low moment of inertia is due to the presence of orbitals at the Fermi 
surface that are very strongly coupled to the deformed potential, which makes the generation of collective angular 
momentum costly. Another case is the smooth terminating bands, which typically have a moment of inertia lower 
than Jrig that is well reproduced by calculations without pairing Q. The reduction is explained by shell effects, i.e., 
gaps in the single-particle level density for certain nucleon numbers. The angular momentum is generated by gradual 
alignment of the individual angular momenta of the particles and holes in incomplete shells, which becomes more and 
more difficult with increasing / causing the decrease in moment of inertia. Without any doubt, there are many more 
examples of substantial deviations of the experimental moments of inertia from the rigid-body value that have been 
reproduced by cranking calculations without pair correlations. We consider these results as further evidence that 
the deviations are caused by the shell structure and not pair correlations. This interesting fact, which has not been 
pointed out enough, will be studied in a systematic way in what follows. 

Fig. I^lshows the experimental ground-state shell energies. They are well reproduced by the calculations in Ref. [l^ 
and other mean-field calculations. The N dependence is governed by the shell structure of the unpaired single-particle 
levels. Both the calculated (Fig. Ill|l and experimental (Figs. 01 and IHl shell moments of inertia show a structure 
similar to that of the shell energies at zero spin, namely: (1) minima at the same magic numbers; (2) higher values 
just below and above the magic numbers than further away; (3) a shoulder or peak around N — 110. The shell 
energies reflect the level density at the Fermi surface. They are negative if the level density is lower than the average 
and positive if it is larger. The relation between the level density and the moment of inertia appears directly in the 
statistical estimate for the moment of inertia (cf. ^) J = g < P >, where g is the level density and < P > the 
average single-particle angular momentum, both taken near the Fermi surface. The Inglis moment of inertia (see 0) 
for collective rotation, X]p/i \^ph\'^ I s-ph, can be estimated by the same expression if we assume that the particle- hole 
energy Cph ^ 1/g and J2ph \^ph\'^ ~< >■ Using for < P > the classical value and for g the Fermi gas value, 
one obtains the rigid-body value |1;]. The deviations from this value are mostly due to the deviations of g from its 
smooth (Fermi gas) value, i.e. due to the shell structure in the level density. This just means it is hard to generate 
angular momentum if the level density is low, so that the moment of inertia is small; and vice versa. Since both the 
ground-state shell energy and the shell moment of inertia are proportional to the level density, their correlation is 
not a surprise. In section HV Gl we shall discuss the relationship between the ground-state shell energies and the shell 
part of the moments of inertia in a more quantitative way. Note that in Figs. IBl 171 and 1121 the shell energy has been 
multiplied by a factor A^^^ where A is the mass of each nucleus. In section HV OI we shall derive this scaling factor, 
which makes the correlation particularly obvious. 

The experimental shell energies at spin 20 (Fig. O look similar to the ones at zero spin. However, the changes with 
N are less rapid. The minima around the main shell closures in {N = 50, 82 and 126) are still present but are less 
pronounced than in the ground state (see Fig. 01. The same is true for the calculated shell energies a,t hu! = 0.3 MeV 
(Fig. I12|) as compared to the calculated ones at zero frequency (not shown) . The reason for the damping of the 
shell structure is the correlation between the zero-spin shell energy and the moment of inertia, which implies an 
anticorrelation between the rotational and the zero-spin shell energy. When the zero-spin shell energy is small the 
rotational energy is large and vice versa. A quantitative estimate will be given in Section If V HI 

Concluding this discussion we claim that the deviations of the moments of inertia from the rigid-body value at 
high spin are determined by the shell structure of a system of independent Fermions confined by a leptodermous 
potential. This implies analogies with non-nuclear systems of confined Fermions. Metallic clusters and quantum 
dots in a magnetic field are examples (see for a review of the analogies), which we shall refer to in the general 
discussion of the shell effects in the next section. The strong deviations of the moments of inertia from the rigid-body 
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value imply that the flow pattern of the mass current must be rather different from the one of a rigidly rotating 
mass distribution. In the case of clusters and quantum dots, the difference between the quantal and the rigid-body 
currents is the electrical current carried by the conduction electrons, which gives rise to an unusually large magnetic 
susceptibility (see e.g., Its name "normal persistent current" alludes to the fact that the large magnetic moment 

is generated by persistent currents like in a superconductor (Mcissner effect), but that the origin of the currents is 
quite different because the system is in the normal state. The analogy is beautiful. The nuclear moments of inertia 
at high spin may be substantially smaller than the rigid-body value as if the nucleus were superfluid. However, it is 
in the normal state and the reduction is caused by the shell effects. In the next section we use the Periodic Orbit 
Theory to obtain more insight into the nature of the shell effects and the mass currents that occur in the rotating 
frame at high angular momenta. 



IV. INTERPRETATION IN TERMS OF THE PERIODIC ORBIT THEORY 



In the preceding section the shell effects in the rotational energy were calculated from the quantal levels in the 
rotating potential using Strutinsky's shell correction method |18| . Now, we want to study them from a different 
perspective. Using the semiclassical Periodic Orbit Theory (POT), which docs not require full quantal calculations, 
we can understand some global characteristics of the rotating many-fermion system. POT has been successfully used 
to explain various aspects of shell structure in fermion systems. The shell structure of a spherical cavity was first 
analyzed by Balian and Bloch 25]. They could relate the spacing between the shells to the length of the shortest 
orbits (the triangle and square) and they predicted a long wavelength modulation which results from the interferences 
between these orbits. This supershell structure was recently found in metal clusters 26]. Strutinsky et al. j27l ] first 
explained how the deformation of nuclei is determined by the classical orbits in the meridian planes of a spheroidal 
cavity. 



A. Basics of the Periodic Orbit Theory 

A detailed presentation of the Periodic Orbit Theory (POT) was given in the book by Brack and Bhaduri 1^. Here 
we review some basic facts needed for our discussion. For a given potential, the level density is decomposed into an 
oscillating part, which represents the shell structure, and a smooth background. The oscillations show up in related 
quantities such as the energy and the moments of inertia, and they are the origin of many structural features of nuclei. 
POT calculates the oscillating part of the level density and of the derived quantities in terms of classical periodic 
orbits in the same potential. More specifically, it aims at the oscillating part of the level density that is averaged 
over a certain energy interval. For a sufficiently wide averaging interval, only the gross structure remains, which is 
described by a few short orbits. This tremendous simplification makes POT a powerful tool for interpreting the gross 
shell structure. We are interested in the gross shell structure of the rotational energy, which has been studied before 
only for the special cases of rotation about the symmetry axis |2^ and the harmonic oscillator [2, 0] ■ 

Let us start with the Fermi gas that describes the smooth behavior without the shell structure. We consider N 
Fermions of one kind in a spherical cavity of radius 

Ro = r,N^'\ (12) 



two of which occupy one orbital state. It is useful to introduce the wave number, k — y^2me/h?, where e is the level 
energy. The relation 

kpRo^l^^y N'^' = bN'/^ 6 = 1.92 (13) 

determines the Fermi wave number kp. Using H12|) . we have kpVo = b. The Fermi energy ep = {hkpY /2m is a 
convenient energy unit. The smooth level density is 

^. . 3Nk 3iVei/2 

5(e) = T7-^ = —T7^- (14) 



2epkp 



The smooth level density in k space is: 



g{k)^dN/dk^g{e){de/dk)^ (15) 

Kp 
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The expressions for the Fermi gas do not depend on the shape of the cavity. They also apply to a deformed cavity with 
the same volume. However, they describe the smooth behavior only correctly in leading order of k. The next order 
(see 113), which accounts for the surface effects, somewhat modifies the relations H13I14(I . These shape dependent 
modifications are not important for the qualitative discussions of the shell structure in the following. 

The level density g is decomposed into the smooth part g and an oscillating part gsh, which contains the shell 
structure, 

9 = 9 + 9sh- (16) 

The oscillating part is given by 

9sh^^gi3 (17) 

where f3 labels the periodic orbits that contribute. 

The shell structure of the levels in a spherical cavity was studied by Balian and Bloch |2^. The classical system 
corresponds to a point mass inside a hollow sphere bouncing elastically from the walls. The periodic orbits are 
polygons in a plane that contains the center of the sphere. Fig. E| (top) shows the simplest, the triangle and the 
square. Each polygon generates a family (3, which consists of all orientations such a planar orbit can have within the 
sphere. The family is characterized by the number of vertices v (number of reflections on the surface) and the winding 
number w (number of turns around the center until the orbit is closed), i.e., /3 — {v,w}. The triangle is {3,1}, the 
square {4,1} and the five-point star {5,2} (see Fig. I14|l . The explicit contribution from each family is given by 

gp{k) = Ap{k) sMLpk + vp)D . (18) 

To be brief, we shall refer to the whole family of such degenerate orbits as the triangle, square, etc. 

As a function of A;, each term in the sum oscillates with the frequency given by the length of the orbital 

Ly^ = 2vsm(py_^Ro, (Pv,w = , (19) 

V 

where if v,w is half the opening angle of one polygon segment. The Maslov index 

1^/3 = l^v^w = ^ + TTW + — (20) 

is a constant phase, which takes into account that each bounce at the surface and each turn around the center changes 
the phase by a constant (see Q)- The amplitude A/i depends on the degeneracy of the periodic orbit: the more 
symmetries a system has, the greater the degeneracy, and the more pronounced are the fluctuations of the level 
density. In the case of the spherical cavity. 



= 2b'/'N'/'^ sin(2^,,^)./^^^^. (21) 
kp * ""^ 

The relative amplitude of shell oscillations compared to the smooth level density ^(e) is ~ (iV)^^/^, that is they are 
comparable for nuclei. The amplitude decreases with the opening angle ify^w of the polygon sections, i.e., it decreases 
as the number of vertices v increases. Hence the triangle and the square have the largest amplitudes. 

The damping factor D{kLi3/jRo) is a decreasing function of its argument. Its concrete form depends on how the 
level density is averaged over k. The wider the averaging interval 7 the more rapidly long orbits are suppressed. In 
the present paper, we do not explicitly average. We consider only the rough dependence of several quantities on the 
particle number. This gross shell structure can be thought of as the result of averaging over a fraction of a shell. 
Hence, it can be understood by considering the interplay of the shortest orbits, which are the triangle and the square 
in the case of the sphere. Note, we have disregarded the diameter orbit, although it is the shortest. It gives only 
a small contribution, because it has a lower degeneracy (two angles are needed to fix a line) than the planar orbits 
(three angles are needed to fix a polygon). Moreover, it does not play a role for the moments of inertia, which are the 
main concern of this paper. 

The shell corrections to the particle number, N^fi, and to the energy, Egi^, are given by the integrals 

Nshikp) = J \sh{k)dk, E,h{kF)^ j ^ ■^{k' -kl)gsh{k)dk. (22) 
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Carrying out the integrations |49j one obtains: 



Nshikp) - - V -Af3{kF) cosiLpkp + vp)D (^^) , (23) 



and 



Esh(kF) = E ( - j ^P^^f) siniLpkp + vp)D j = E \-) 9p{kF). (24) 



The period of revolution of a particle moving with the Fermi momentum hkp on the orbit [3 is given by 



h_ ^ fekp_ ^ ^epRo^.i/s^ (25) 
Tf3 mLp b Li3 

The shell energy is Egh ^ epN^^^ as compared with the smooth part of the Fermi gas, which is ^ epN . 

Eqs. H23I24|) should be understood in the following way. Both Ngh and Esh depend on the particle number via 
Ro = ToN-^/^, where we have introduced the "smooth" particle number for a clean notation. The actual particle 
number N = N + Nsh{N) and the shell energy Esh{N) are functions of the parameter TV, which then defines Esh{N) 
in a parametric form. However when discussing the extrema of Egh, we need not resort to this sophistication, because 
C0B{Lpkp + vp) ~ where s\n{LpkF + V13) = ±1- 



B. Basic shell and supershell structure 

Let us discuss the main features of the shell structure of the spherical cavity 'g^ as an educational example. The 
most important orbits are the triangle (/3 — A) and the square (/? = □), to which we restrict the sums over /3. They 
are the shortest orbits with the lengths La — 5.19i?o and ~ b.QQRo- Since, La ~ Ln, one has ta ~ « r, 
Aa ~ ~ A and _Da ~ -Dn ~ D. Using the addition theorem for the sine function one finds 

EsH = 2(^^^ Asm{kFL + D)cos{kFAL + Aiy)D, (26) 

with L = (La + Lq)/2 = 5.42_R and AL = (Lq — La)/2 = 0.24i? and the analogous definitions for D and Ai'. A 
well known phenomenon is encountered: the superposition of two oscillations with similar frequency results in a beat 
mode. The fast oscillation represents basic shell structure and the slow beat mode was called supershell structure 

M- _ 

The phenomenon of supershell structure was observed in Na clusters |2g , almost three decades after it was predicted 
psf. It is realistic to assume that the conduction electrons move in a cavity. The basic shell closures correspond to 
the minima of the sine function, which lie at kpL + P ^ 7r(2n + 3/4). With L = 5.42i?o and kpTo = b one finds 
that the closed shells appear at N^^^ = 0.60n + c. Overall, the experimental magic numbers in Na clusters are well 
reproduced by the relation N^/^ = 0.61n + 0.50 The magic numbers 58, 92, 138, 198 give 921/^ - 58^/^ = 0.64, 
138^/'^ — 92^/'^ = 0.65, and 198^/"^ — 138^/'^ = 0.66. The somewhat longer period for the small clusters as compared 
to the POT value for the cavity is due to the assumption of a cavity instead of a potential of finite depth and surface 
thickness. 

The slow oscillation is the supershell structure and has a half period of kpAL = tt, which corresponds to L/2AL « 12 
shells. In the experiment, the beat minimum appears around n = 15 26J. When the cosine function of the slow 
oscillation changes sign, the maxima of the fast oscillation become minima, i.e., the new shell closures, which are 
shifted by half a shell as compared with the ones in the lower beat. This phase shift was also observed at n = 15 [26|. 
A more careful application of POT than given here (see ^,2^ 3Q]) removes the discrepancy between the calculated and 
observed shell number where the beat minimum appears. 



C. Deformation 



In the middle between the closed shells, Esh > for spherical shape. Nuclei and alkali clusters reduce this shell 
energy by taking a non-spherical shape, i.e., they avoid a high level density near the Fermi level. This is analogous to 
the Jahn- Teller effect in molecules. Due to symmetry, the electronic levels of a molecule may be degenerate. If such 
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a degenerate level is incompletely filled, the molecule changes its shape such that the degeneracy is lifted. Whereas 
the final shape of the molecule is determined by the balance between this driving force and the restoring force of 
the chemical bonds, in the case of nuclei and clusters a shape is attained that minimizes the level density near the 
Fermi surface. The corresponding gaps in the single-particle spectrum are referred to as " deformed shells" lU • The 
optimal shapes are described by the few lowest multipoles. This has been known for a long time for nuclei, where it is 
experimentally confirmed. Also in the case of alkali clusters, several multipoles are needed to describe the equilibrium 
shapes (see, for example the calculations in [s^). The calculations of the shapes in the present paper include the 
multipoles necessary for a completely relaxed axial shape. POT permits us to understand some basic features of 
the deformation. However, a full understanding of the interplay between the multipoles has not yet been reached 
(concerning octupoles see p3^). 

Our present interpretation of deformed shapes in terms of POT is based on the analysis of the spheroidal cavity. 
Strutinsky and coworkers did the pioneering work |27| by discussing Egh as function of the particle number N and 
the ratio rj of the long and short semi-axes, which are, respectively: 

Ri = roN^'^rf/^ R, = roN^'^rj-^'^ prolate (27) 

and 

Ri = ToN^'^r^^''^ Rs = roN^I^'q-'^'^ oblate. (28) 

The volume of the cavity is deformation independent. Fig. 1151 (tor)) shows the shell energy contours Esh as function 
of the deformation parameter a (see H^), calculated for the Woods-Saxon potential with a very thin diffuseness 
(a=0.05) and no spin-orbit coupling, which is practically identical with the cavity. In this case, to reach the same 
accuracy the integration mesh for calculating the Hamiltonian matrix elements was taken w 10 times finer than for 
the usual calculation with the parameters from the Tab. For the considered deformation range the two deformation 
parameters are related by 77 « \/(l -I- |a|)/(l — 

The most important orbits are the polygons lying in the meridian planes, because they have a two- fold degeneracy: 
all polygons in one meridian plane have the same length and all meridian planes are equivalent. Again, the triangle 
and the rhombus are the most important orbits, which are shown in Fig. 1131 As for the sphere, their interference 
generates a beat pattern [3 US 7 i-e., the shell correction E^h has the form (|21l, where A is difi^erent. The basic 
shell structure is governed by sui{kFL± + v), where L±_ = (La + L^)/2 and _L indicates the meridian plane. (The 
reason for the labeling becomes clear below in the context of rotation.) In order to keep the expressions simple, we 
approximate Lj_ « L^, following Ref. [23 ■ The length of the rhombi is: 

The equilibrium shape corresponds to the minimum of Egh, i-e. to 

3 

kpL^ + u ^TT{2n+ -). (30) 

These lines, L±{N,ri) = const, are also shown in Fig. El As seen, the valleys and ridges follow the lines, L±{N,r]) = 
const. They are nearly horizontal on the oblate side, because the function yrf^ is roughly constant in the 

interesting range of rj. However, on the prolate side, constant length corresponds to approximately N cx ?7~^, which 
results in the down-sloping curves. If one starts from a closed shell and spherical shape [j] = 1 or a — 0) taking 
particles away, it is energetically favorable to follow the valley on the prolate side. This is H. Frisk's explanation |3J| 
for the preponderance of prolate over oblate nuclei. The smooth increase of the deformation with decreasing N is 
another experimental fact (see e.g., ISlj ) which is explained by the down-sloping of lines, L{N,r]) = const, on the 
prolate side. If N decreases further, the spherical shape in the next lower valley eventually becomes energetically 
favored. Since the two valleys are separated by a ridge, the deformation decreases abruptly, when the valley near 
spherical shape takes a lower energy. This explains the sudden onset of nuclear deformation when the open shell is 
entered. Both features have been first pointed out in Ref. . The traverse of the meridian ridge is clearly seen in 
the shell energies of Na clusters calculated in Ref. [s^l (cf. Fig. 2 therein, minimization with respect to a only). 

Let us elaborate on the above analysis of deformation from previous work. The N ~ rj landscape contains more 
structure than the sequence of ridges and valleys generated by the meridian orbits. It shows an interference pattern 
with a second set of valleys and ridges that is generated by the orbits that lie in the equator plane and its neighborhood. 
The equator orbits are regular polygons as in the spherical cavity. In contrast to the meridian orbits, they are only one- 
fold degenerated with respect to a rotation around the center. Therefore the valley-ridge structure is less pronounced 
than for the two-fold degenerated meridian orbits. Again, the valleys and ridges correspond to a constant length 
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L|| = (La + Lu)/2, where \\ indicates the equator plane. For simphcity we approximate _L|| « L\j, which on the 
prolate side is given by 

Lu = 4%/2roiVi/\-i/^ (31) 

Fig. Elalso shows the lines, L^\{N,'ri) — const. The ridges start at the maxima of the spherical cavity and the valleys 
(not shown) at the respective minima. This is analogous to the meridian orbits, however the slope of the lines of 
constant length is positive, since N oc rj. Both types of orbits originate from the same family at spherical shape. The 
equator ridges and valleys modulate the bottom of the meridian valleys. If one starts from a closed shell at 77 = 1 
taking particles away and moves along the bottom of a meridian valley one passes a hump that is generated by an 
equator ridge and continues into a depression generated by an equator valley. The depression is the deformed shell 
closure. There, the meridian and equator orbits interfere constructively. The equator hump and depression are clearly 
seen in the shell energies of Na clusters calculated in Ref. j32| (cf. Fig. 2 therein, minimization with respect to a 
only). 

D. Superdeformation 

Matsuyanagi and coworkers [3^ studied the relation between the periodic orbits and the appearance of a pronounced 
shell structure in a prolate cavity with an axes ratio of about 2:1, which corresponds to the nuclear superdeformed 
shape. They found that different orbits are responsible for the shell gaps causing superdeformation. The meridian 
orbits lose importance with increasing deformation, because they become longer. The radius of the equator plane 
shrinks with increasing deformation. Therefore, the equator orbits become shorter and more important. Most 
prominent are the orbits that close only after two turns around the symmetry axis (w = 2 in Eq. H19|) ). The 
five-point star is the simplest orbit of this type, which is shown in Fig. ^1 Between 77 = 1.5 and 2, these planar orbits 
become unstable. They give birth (bifurcation) to three-dimensional orbits with the following path: during the first 
turn the orbit is above the equator plane and during the second turn it is below. Pictures of these three-dimensional 
orbits can be found in Ref. [33 • The five-point star and the double-traverse triangle determine the superdeformed 
shell structure. The butterfly orbit shown in Fig. 1141 is prominent in the deformation interval 1.4 < rj < 2. It bifurcates 
from the double-traverse diameter orbit in the equator plane. These orbits become more important than the shorter 
orbits with only one winding around the center, because the amplitude strongly increases near the bifurcation points. 
In fact, the amplitude goes to infinity if one uses the standard stationary phase approximation. A quantitative 
description of the enhancement must be based on a more accurate treatment, which Magner et al. worked out for the 
spheroidal cavity [37l| . 

The lines of constant length of the five-point star orbit in the equator plane are also included in Fig. The change 
from the low- deformation to the high-deformation shell structure is more clearly seen in figures that extend to larger 
deformation (see Refs. '27','3^). However, some of the superdeformed shell structure is visible on the right fringe of 
the figure. The fact that double-traverse orbits are responsible for the superdeformed shell structure, is most directly 
reflected by the shorter wavelength of the oscillations of the shell energy. The period of the spherical shell structure 
is ^ 27r/(5.4 b) = 0.60 (cf. Sect. HV^ . The double triangle has the length 10.4 roA^i/^^-i/s^ which gives a 

period of AN^^^ = 27r77-'^/^/(10.4 b) = 0.40 for rj — 2. In accordance with this estimate, the quantum calculations in 
Refs. [23,13^ give a ratio of about 1.5 for the periods of the spherical and superdeformed shells. 

A closer inspection of Fig. shows that at normal deformation (a « 0.3) the ridges and valleys from the meridian 
orbits seem to be modulated by both the single- (triangle, square) and double- (flve-point star, double traverse triangle) 
traverse orbits in the equator plane. Thus, our discussion in section llV Cl oversimplified the interplay between meridian 
and equator orbits by disregarding the double-traverse orbits, the inclusion of which accounts for details of the N — rj 
landscape in the upper panel of Fig. (e. g. the bump at = 90 and a — 0.3). 

E. Application to nuclei 

The cavity model differs from nuclei in three important aspects: i) the nuclear potential does not jump to infinity 
at the surface; ii) there is spin-orbit coupling; and iii) there are protons and neutrons, which fill the shells differently. 
We will discuss each of these briefiy. 

i) Refs. 27, 30] compared the shell energies of the cavity with the ones of the Woods-Saxon potential. They found 
that the qualitative features agree rather well. We are using as a cavity-like model the Woods-Saxon potential with 
the same depth, but with very small diffuseness (a = 0.05) and no spin-orbit coupling. The results for such a model 
are close to the true cavity model. 
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ii) Incorporating the spin-orbit coupling into POT (cf. [s^ Is^ lioj ') is complicated because the spin cannot be 
treated semiclassically. Simple interpretations as for the spinless POT have not been found yet. However, the shell 
energies calculated from Woods-Saxon potentials with and without a spin-orbit part show strong similarities if one 
compares them at the same fractional filling of the spherical shells, as seen by comparing the upper and lower panel 
of Fig. El In the following we will assume that this also holds for the rotational energies. 

iii) Depending on the nuclide, protons and neutrons may fill a different fraction of the shell and contribute in a 
different way to the net shell structure. We show the experimental data and the calculations as functions of the neutron 
number because the neutron shell structure is more clearly visible. Around the spherical shell closure at A'^ = 126, 
Z is not too different from 82, i.e., the proton and neutron shell contributions are in phase. They get progressively 
out of phase with decreasing N. At N = 82, Z is around 64, i.e., about mid shell. For 50 < iV < 82 the protons 
are out of phase with the neutrons. The neutron shell contributions are stronger than the proton contributions in 
the heavy nuclei. Therefore, the neutron shell structure shows up in the N dependence of the energies and moments 
of inertia. However, it is somewhat modified due to the proton shell contribution, which depends on Z. For our 
qualitative discussion, we may compare the total shell contributions in nuclei at a given fractional filling N of the 
neutron shell with the ones calculated for the cavity at the same fractional filling TV. For N > 100, where the neutron 
shell contribution are not too much out of phase with the proton contribution, one expects that the N dependence of 
the shell effects is about the same as in the cavity that contains only one kind of fermions. For N < 100 one expects 
a weaker N dependence, because the protons no longer enhance the neutron shell structure. In fact, the proton 
shell structure may counteract and modify the neutron shell structure by changing the deformation. A detailed 
discussion of this interplay goes beyond the scope of this paper. However, some features of the shell structure that 
are a consequence of this interplay will be discussed below. 

The magic numbers 50, 82, 126 (and 184 predicted) correspond to periods of the spherical shells 82^/'^ — 50^^'^ = 0.66, 
126^/'^ — 82^/-^ = 0.67, 184^/'^ — 126^/'^ = 0.67, which are the same as in the small clusters and somewhat longer than 
the period AA^i/3 = 0.60 given by POT for the cavity (cf. section HV B|l . The closed superdeformed shells appear 
at A^ = 112 and 88, which corresponds to a period of 112^/'^ — SS^/'^ = 0.37. The POT period of the superdeformed 
shells in the cavity is AN^/^ = 0.40 (cf. IIVDI) . The ratio of 1.7 between the experimental periods of the spherical 
and superdeformed shells in nuclei is somewhat larger than the POT prediction of 1.5. The ratio reflects the different 
lengths of single-traversed orbits in the sphere and the double-traversed orbits in the equator plane of the spheroid 
with the axes ratio 2:1. 

The experimental shell energies at zero spin in Fig. Elshow quite clearly the interference between the equator and 
meridian orbits at normal deformation. Going down from the closed shell at A = 126 one climbs up the bottom 
of the valley generated by the meridian orbits, crossing the ridge generated by the equator orbits, and reaches the 
deformation region generated by the constructive interference of both the meridian and equator valleys (N=98). One 
then follows the valley generated by the equator orbits. Along the bottom of this valley one has to go over the 
ridge generated by the meridian orbits (N = 90) in order to reach the spherical minimum at A^ = 82. Though less 
pronounced, the same pattern is seen in the shell 50 < A^ < 82. 

The Fermi gas relations are modified for the nucleus because each orbital state is occupied by four fermions 
(proton, neutron, spin up, spin down). The radius is i? = TqA^^^ with Tq = 1.2 fm. Assuming N = Z = A/2 we have 
kpTo — and ep — 34 MeV. Hence, for the action one has 

/4 \ 1/3 

kFRo = kFroA^^'' ^bi-j ^hN^I^, (32) 

as in the case of a cavity with only one kind of fermions. Studying the rotation below, we encounter the inverse energy 

2mif ^ 2m {rpA^/^kpf _ {h{A/2Y'^f _ {bN^/^f ^^^^ 

tfi ~ [nkpy 

In the heavy nuclei, where A^ > A/2, there appears the ambiguity, whether one should use for N the neutron number 
or A/ 2. The difference in the estimates is not significant in the context of our qualitative discussion. In order to be 
definite, we use A^ = A/2 whenever we refer to A and the actual neutron number whenever we refer to TV. 

F. Influence of rotation 



Rotation is taken into account by applying POT to the Routhian 

H' - w/, 



(34) 
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where H is the Hamihonian of the deformed cavity, lo the angular velocity and I the projection of the orbital angular 
momentum on the axis of rotation. The spin is disregarded. If 00 is interpreted as the Larmor frequency, the Routhian 
H34|l agrees with the Hamiltonian of a system of electrons in a weak magnetic field. Hence, we can directly use several 
results from studies of a cavity in a magnetic field. 

The smooth part of the moment of inertia is the rigid-body value [soj , which is for the spherical cavity containing 
only one kind of fcrmions 



2t2 



Jr.n = -mR'N = 



In the case of nuclei it becomes 



Jrtg = -mR^A = 



22/35, 



.^5/3^ 



(35) 



(36) 



The classical orbits in the cavity are modified by the rotation. The particle moves on a curved trajectory between 
the reflections on the surface. At the Fermi level, the deviation from the straight line is proportional to the ratio of 
the velocity of the cavity Vr = tlo and the particle velocity in the non-rotating cavity vp — PF/Tn. Using the maximal 
value of Vr at the surface, we have 



VF 



mRoo 

PF 



mToN^^^LO bfko 



PF 



2eF 



(37) 



where we used the Fermi gas estimates. From the maximal angular momentum Imax we consider in this paper (c.f. Fig. 
IHJ and the experimental moments of inertia J^exp one has fiio < Imax/ Jexp, which means that the ratio vr/vf < 0.1. 
Therefore, the Cranking term col can be treated in a perturbative way. Kolomietz et al. applied the perturbation 
theory to a rotating spherical cavity. Tanaka et al. E3| and Frauendorf et al. |2^ used it for studies of the magnetic 
response of electrons in a spherical cavity. Creagh pa | formulated the perturbative approach in a general way. The 
method and the treatment of a weak magnetic field are exposed in detail in Q , which we follow here. 
The change of the action due to rotation is given in first order by 



ASb = 



) / Idt — UJTplp = h^i3. 



(38) 



The integration runs over the unperturbed orbit in the non rotating cavity. In the case of the spherical cavity, the 
angular momentum Ip of the orbit is conserved and the integration is trivial. In case of the spheroidal cavity, is 
conserved for the equator orbits but not for the meridian. For the latter, Ifj is the average angular momentum of the 
orbit which must be found by evaluating the integral (|38|) . Since loI — •m{r x v) ■ lo, 



(f X ds) ■ uo — 2m / df 

13 Jl3 



2mApLi! cos 6 = h^p cos 9, 



(39) 



where A/^ is the area enclosed by the orbit and 6 is the angle between the normal of its plane and the axis of rotation. 
Hence, $ is the "rotational flux" in units of h, of the vector field 2mL0 enclosed by the orbit. In the analogous case of 
particles in a magnetic field, $^ is the magnetic flux in units of the elementary flux quantum. 

The rotation manifests itself in the appearance of an additional modulation factor M in the expressions for the 
level density (|18|) and the shell energy H22|) . The modulation factor is given by the avera ge o f exp(i<I>) over all orbits 
of the same length, which belong to the family (3. In the case of a spherical cavity l42| . the modulation factor 
becomes 



(40) 



where $^ is the flux through the orbit perpendicular to the rotational axis (or the magnetic field) and jo the spherical 
Bessel function. 

In the case of the spheroidal cavity, we distinguish between the rotational axis being parallel or perpendicular to 
the symmetry axis. For parallel rotation, only the equator family carries rotational flux. It consists of all polygons 
rotated around the center. Since $ is independent of the rotational angle one has 



A^||($/3) = cos(<I>;5j = cos{2mApuj/h) = cos{ujTplp/h). 



(41) 



The cosine function appears because in each family the orbits appear in time reversed pairs (the particles run clockwise 
and counterclockwise on the orbit). The case is analogous to spherical quantum dots (two-dimensional circular 
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potential pockets that confine electrons on a surface) in a perpendicular magnetic field. The modulation factor has 
been observed as oscillations of the electric current through the dot, which oscillates with the magnetic field strength, 
where the period of the oscillations is given by the area of the triangle (see [sl li^). 
For the meridian orbits 



MA'fp) = ^ r e'^f'^°^'de = Jo{<^p), (42) 

where Jo is the Bessel function and the flux through the orbit perpendicular to the axis of rotation. Since the area 
enclosed by the family members lying in one plane differs for the meridian orbits one must average the modulation 
factor over all these orbits (see below). Fig. 1161 shows the three types of modulation factors. 

In the analysis of the experiment and the microscopic calculations we assume that the rotational energy has the 
form which corresponds to an expansion of Ai up to second order in ^p. Using the expansions of cos(<i>^), 

Jo($/3) and joi'^p), one finds 

M^l- a*^, (43) 

where, respectively, a||=l/2 and a^=l/4 for the equator and meridian orbits in the spheroidal cavity and qq—I/Q for 
the spherical cavity. One may derive this result directly by expanding exp(i$^cos0) before integration. The leading 
term is quadratic, because the linear term becomes zero when averaging. The next term is (<1>^ cos0)'*/24, which can 

/ N 1/4 

be neglected in our qualitative discussion as long as '^p < 1.5/ I (cosO)^ j . As seen in Fig. ^1 1 — A4p deviates from 

its quadratic approximation by less than 35% for <i>^ < 2, where the quadratic approximation is better for Jo($^) 
and joi^p) than for cos(<l>^). 

Inserting the quadratic approximation H43|> of the modulation factor into eq. (|24|l one obtains: 



E'^^ikp) = E^nikp) ~ -jJshikFl (44) 



and 



(45) 



This low-o; version has been used in by Frauendorf et al. |2J| for analyzing the magnetic response of spherical metal 
clusters and circular quantum dots. In the present paper we have added the meridian orbits. 
The angular momentum of the orbit is 

^ 2mAp ^ 2hkFAp ^ 2ft^^^i/3^ ^^g^ 
T/3 RLfj 

where Ap — \/~A^ and A^ denotes the average over the degenerate orbits in one and the same plane. The exact form 
of this average will not be derived in this paper, because it is not important for our qualitative discussion, as seen 
below. Working it out, would demand a substantial sophistication of POT. 
The area of the regular polygons in one of the planes of the sphere is 

= A,,^ = - sin )B? = - sin r^TV^/s. 47 

2 V 2 V 

The area of the polygons in the equator plane of the prolate spheroid is given by the same expression containing 
instead of R. Averaging over the different orientations in the equator plane does not give anything new because the 
area is constant. For the basic shell structure, we take the mean of the contribution from the triangle and the square, 

,43, 

where || indicates that the rotational axis is parallel to the symmetry axis. In the case of the meridian orbits, the 
areas are rather complicated expressions in terms of the half-axes (see 0|). Let us consider the rhomboidal orbits. 
The orbit, whose diagonals coincide with the axes of the ellipse, encloses the area 



","3:*: 



(49) 
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The area enclosed by the orbit consisting of Hnes parallel to the axes of the ellipse differs from (|49|l by the factor 
2/\/l + '7 + l/'7- The areas of the other rhombi fall in between these limits. Since for 77 < 1.5 the differences between 
the areas are less than 2%, we use ■ We assume that the area of the triangles changes with deformation as given 
by Eq. (|49|l and use the approximation 



^1/3 prolate ^^^^ 
7^-1/3 oblate ' 



where _L indicates that the axis of rotation is perpendicular to the symmetry axis. 

Let us now discuss the accuracy of the quadratic approximation for the modulation factor in the case of nuclei, 
which is illustrated in Fig. ^| The rotational flux in the sphere is (cf. (|33I39I48|| ) 

^ 2m(lMRl)hw 1.0362 

.2 = ■ 51 

ep 

In the microscopic calculations, we use hoj = 0.3 MeV, which corresponds to fluxes of $ = 0.6 and 1.1 for A ~ 80 and 
200, respectively. Hence, the quadratic approximation is rather accurate. In the case of the experiment, we extract 
the shell contributions to the moment of inertia mostly from the last lOh, i.e., the average spin / « Imax — 5. Using the 
experimental moments of inertia, we may estimate the frequency as huj ~ I / J^exp, which gives hio — 1.0, 0.5, 0.4 MeV 
and <I> = 2.0, 1.6, 1.5 for the regions A = 80, 160, 200, respectively. Hence, for the bulk of the data the quadratic 
approximation is not very accurate, but still acceptable for our first systematic analysis. It becomes problematic for 
the light nuclei and for some the cases that reach very high spin (cf. Fig|Sl). Including the deformation reduces the 
flux for rotation about the symmetry axis. With 77 = 1.3, one finds (f>|| = 1.34 for A = 160, for which the quadratic 
approximation of cos(<l>) becomes quite reasonable. The deformation increases the flux for rotation perpendicular to 
the symmetry axis. With 77 = 1.3, one finds — 1.8 for A — 160. A parabola that approximates A^($) at large 
values of $ has a lower curvature than the parabola that approximates the low-$ part (see Fig. 116(1 . Hence, the 
deviations of Ai{'^) from the quadratic form will tend to reduce the shell contributions to the moments of inertia that 
we derive from the data. 



G. Moments of inertia 

The moments of inertia are given by: 



and 



Jsh\\ = llAikp) &in{L\\kF + v)D (^^^ , (52) 



Jshi^ = ^liAikp) siniL^kp + D)D ( ) . (53) 



The ratio (cf. Eqs. (|i^ll46l^6|l ) 

^ ^ l^AepN-^/' ^ (^1/3)2^5/6^-5/3 ^ ^-1/6 (54) 



is of the order of one, i.e. the shell contribution to the moment of inertia is comparable with the smooth rigid-body 
value. 

One may rewrite Eqs. H52l53fl as follows: 

Jsh\\ = ijashwikp) (55) 

and 

Jsh± = ^l'^^gsh±{kp). (56) 

This is a more quantitative statement of the relation between the level density and the moment of inertia, which we 
discussed in Sect. IIII Bl The deviation of the moment of inertia from the rigid-body value is not determined by the 
total deviation of the level density from the Fermi gas value but only by the part generated by the orbits carrying 
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rotational flux. This is indicated by the subscripts \\ and _L for the parallel and perpendicular orientation of the 
rotational axis. 

Taking into account Eq. H24|l . one may also relate H52I53(I to the ground-state shell energy: 

JsMI = (57) 

and 

Jshi^ = ^{klA^fEsh±, (58) 

which gives Jsh ~ {^)^ N^'^'Esh ~ 400 tiev- N^'^'Esh « lomiiev- ^^'^'^sh- This is why in Figs. H [7| and [El 

the shell energy was scaled by a factor A^/'^. The scaling factor accounts for respective scales on the plots of Jsh and 
Esh, which is remarkable considering the qualitative nature of our arguments. As in the case of the level density, only 
the part of the shell energy originating from the orbits enclosing rotational flux contributes to the shell moment of 
inertia . 

The terms in the sum (145(1 giving Jsh contain the additional factor {TplpY (as compared to the terms in the sum 
(|24|1 giving Egh) which is proportional to A^, the square of the area enclosed by the classical orbits (cf. discussion in 
the preceding paragraph). This factor tends to favor the contribution of the longer orbits to the moment of inertia as 
compared to the shell energy, which will introduce more fluctuations in the variation of the shell moment of inertia 
as a function of neutron number. These fluctuations would be damped by averaging over nuclei (which we have not 
done) because the damping factor in the sums would suppress the longer orbits. 

The relations (|57I58|I allow us to discuss the gross shell structure of the moments of inertia with the help of Fig. 
[T51 The shell structure of the moment of inertia was first discussed in 0. Fig. 2 of Ref. j3] shows the neutron part of 
J7j_ for a = 0.3. There is a maximum below N = 90, a minimum near N = 106, and another maximum at = 120, 
which correspond, respectively, to the ridge, valley, and ridge, generated by the meridian orbits in Fig 1151 fbottom). 
Fig. 3 of the same paper shows the proton part of J± for a = 0.3. There is a maximum at Z = 58 a minimum at 
Z = 68 and another maximum above Z = 80, which correlate well with the sequence ridge, valley, ridge, generated 
by the meridian orbits in Fig. 1151 fbottom). The interpretation of Fig. II II and its experimental counterpart Fig. ^is 
more complicated, because the deformation and the orientation of the rotational axis are optimized. For this reason, 
we first discuss the four subsets of calculations shown in Fig. ^| in each of which the axis of rotation and the sign of 
a are fixed, with the help of Fig. El 

For prolate deformation and rotation perpendicular to the symmetry axis the meridian orbits determine the shell 
moment of inertia, because only they enclose rotational flux. Starting at the iV = 82 minimum and following the 
east path of lowest elevation around the N = 100 mountain, one goes over the meridian pass, which shows up as 
the maximum of the Jsh±, arriving at the deformed minimum around N ~ 98. Then one follows the valley to the 
N = 126 minimum, realizing that there is no TV = 106 ridge, which is generated by the equator orbits. Accordingly, 
Jsh± decreases until the minimum at A'^ = 126 is reached. 

For rotation about the symmetry axis, Jsh\\ is determined by the equator orbits. Following the same path, one goes 
through a valley to a minimum at A^ = 98 and over a hump at A^ = 106 (the ridges generated by the equator orbits are 
dashed in Fig. which determine the minimum and maximum of Jsh\\, respectively. The behaviour near N = 82 
will be explained below. For oblate deformation one goes on the west path around the A^ = 100 mountain (Fig. I15|l . 
where a reaches values between -0.15 and -0.20. The maximum of Jsh± at A^ ~ 105 is caused by the meridian ridge 
and the maximum of Jsh\\ sA. N — 90 reflects the equator ridge. 

All four calculations look similar around N = 82. There is a minimum at A" = 82 enclosed by two spikes at A^ = 80 
and A^ = 86, 88. Inside this interval the shape is spherical. When A" goes outside the interval, the deformation 
suddenly jumps to a substantial value and the shell contribution drops. The deformation jump is a consequence of the 
topology of the map, i.e.. Fig. E| Near the closed shell, the equipotential lines are ellipse- like with the center at the 
spherical minimum. For given A^, the energy minimum lies at a — 0. In the middle of the shell, the equipotential lines 
become ellipse-like around the mountain top. Hence, when A^ increases from 82, a stays zero as long as the contour 
lines curve downwards. When the curvature changes sign, a jumps to a substantial value {\a\ ~ 0.15). The analog 
happens when A^ decreases from its magic value, 82. Hence within the interval 80 < Af < 88, the moments of inertia 
reflect the parabolic A^-dependence of the shell energy for spherical shape. Outside, the deformed shell structure shows 
up. The distinct structure around A?^ = 82 is due to the interplay between the proton and neutron shell contributions, 
which are out of phase. The isotones around N — 82 have 60 < Z < 70, i.e., they have a large positive proton 
shell energy for a = 0. Accordingly, the protons make a large positive shell contribution to the moment of inertia, 
which mostly compensates the negative neutron contribution. Therefore, the total shell contribution to the moment 
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of inertia is only slightly negative at iV = 82 and becomes positive on both sides. When the deformation jumps to 
finite values, the positive proton shell contribution is drastically reduced, which ends the spikes. 

Around N = 126, the protons and neutrons are in phase and the total shell contribution stays negative. The 
moments of inertia reflect the shell energy. Going down on the prolate side, the ridge due to the equator orbits shows 
up in J^shW the bump near N = 106 in Fig. Uni which is not seen in J^sh± , because the equatorial orbits enclose no 
flux. Going down on the oblate side, the maximum of J7s/i± at N ^ 100, reflects the upper side of the mountain in 
the shell energy, which is generated by the meridian orbits. The maximum of JTsh\\ at w 90 reflects the lower side 
of the mountain in the shell energy, which is generated by the equatorial orbits. 

For N increasing from 126, the calculations assuming a > result in a spherical equilibrium shape. The same 
is true for the calculation constrained to a < and perpendicular orientation of the rotational axis. In these cases 
the shell moment of inertia has the deep minimum at TV = 126 characteristic for spherical shape. In contrast to the 
N — 82 region, where protons and neutrons are not in phase, the proton number here is Z 82, in phase with the 
neutrons. The combined contribution to jTsh. leads to the strongly negative values, which reach almost the rigid-body 
value {Jcaic becomes nearly zero). The calculation that assumes a < and parallel orientation of the rotational axis 
gives —0.15 < a < —0.20 for most nuclei. The shell moments of inertia are much less negative than for a spherical 
shape, as the comparison with the other calculations in Fig. 1101 shows. Only for Z = 82 and iV < 120 is the shape 
spherical, and then the shell moment of inertia decreases strongly. We discuss the competition between the shapes 
and different orientations of the rotational axes in section llV HI 

Fig. ^2 is a combination of Figs. ^| Near TV — 82, it shows the two spikes at = 78 and A^ ~ 86, which are 
caused by the sudden transition from spherical to deformed shape. For 90 < A^ < 100 the shape is prolate and the 
rotational axis perpendicular to the symmetry axis. The shell moment of inertia Jsh^ gradually decreases from about 
20 h^/MeV to -10 T?lMeY. Above N = 102, rotation about the symmetry axis is preferred, and the peak of Jsh\\ 
A^ = 106 shows up. Above A^ = 110, the shape is oblate and the rotational axis parallel to the symmetry axis. Only 
the Z = 82 chain has a spherical shape aX huj = 0.3 MeV. 

The experimental shell moments of inertia have a similar A^ dependence as the calculated ones. The A^ = 106 
peak is less pronounced than in the calculations. This can be understood as follows. In the calculation it is assumed 
that the rotational axis is parallel to the symmetry axis, i.e., the yrast line is a sequence of high-K band heads. The 
experimental yrast lines contain only few such states. Most of the states entering the fit arc band members with 
I > K. For such states, the rotational axis is neither parallel nor perpendicular to the symmetry axis, rather it 
has an intermediate orientation. As a consequence, the shell moment of inertia should be between J^sh\\ ^-nd J^shJ., 
i.e. smaller than J^sh\\- In our analysis of the experiment, we subtract the rigid-body moment of inertia for rotation 
perpendicular to the symmetry axis. For a tilted rotational axis the rigid-body value will be smaller and jTs/i larger. 

H. Shell energies at finite spin 

Fig. El shows the shell energies of a cavity rotating perpendicular to its symmetry axis. Lets us first consider 
spherical shape {a = 0). The modulation factor A^q($) given by eq. (|40|l damps the shell structure of the non- 
rotating cavity. The flux through the orbits that govern the basic shell structure (mean of the triangle and square) 
is 

$ = IMb'^N'^/^hio/eF = 0.178 {MeVy'^N'^^^huj, (59) 

which gives 1.26 for A^ = 115 and huj — 0.3 MeV. As seen in Fig. ^1 the modulation factor A^q(1.25) — 0.76. The 
top of the A^ = 115 mountain is at 18 MeV for io = (cf. Fig. I15|l and 11 M eV for hui = 0.3 MeV, which corresponds 
to a reduction of 0.6. For huj = 0.6 MeV one has $ = 2.5 which gives A^q(2.5) = 0.24. The above maximum energy 
is 2 MeV for this frequency, corresponding to a reduction factor of 0.1. Similar qualitative agreement holds for the 
other minima and maxima at a = 0. In the case of the realistic potential (cf. Figs. 1151 and [THll . the reduction of the 
shell energy in the quantal calculation is 0.81, 0.70, and 0.73 for huj — 0.3 MeV and A^=82, 100, and 126, which is to 
be compared with the POT estimates 0.84, 0.79, and 0.73, respectively. 

Comparing the experimental shell energy at / = 20 (Fig. [TJ with the ground-state shell energy (Fig. O, one sees 
that the depth of the minimum at neutron number A^ — 126 is reduced by about a factor of 2. A flux of $ = 1.9 is 
needed to make the modulation factor equal to 0.5. The estimate of the rotational frequency is not very accurate, 
because of the irregular sequence of the yrast levels. Near N = 126, the experimental moments of inertia scatter 
between 30 and 45 fi^/MeV, which gives ha; = 0.66 - 0.44 MeV and $ = 2.95 - 1.96 for / = 20. Another possibility 
to estimate the frequency is to average the level spacings {E{I) — E{I — 2))/2 over several spins near I — 20. The 
sample of nuclides, for which enough data exist, is (Z, A^) = (82,120), (82,122), (86,123), (86,126), (86,128), (88,126), 
(88,128), (88,129), and (88,130). The distribution of the frequencies is fuv = 0.36 ± 0.22 MeV, which corresponds to 
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the flux $ = 1.6 ± 0.9. Hence, the estimated frequencies and reduction of the shell structure are consistent. However, 
more data are needed to improve the statistical significance. 

Let us now consider a substantial deformation. Only the orbits that carry flux enter the modulation factor, and 
their contribution to the shell energy changes with to. The contribution of the other orbits remains the same. This 
simple observation is the key for understanding the evolution of shell energies with frequency. Figs. 1 1 71 and [T51 flower 
panels) show the shell energy for rotation perpendicular to the symmetry axis. With increasing frequency lu, the 
valley-ridge system generated by the meridian orbits is attenuated and mostly gone at hio = 0.6 M eV. What remains 
is the upsloping valley-ridge system generated by the equator orbits. For — 100, a = 0.3, and fku = 0.6 MeV, the 
flux ^± = 2.5 (cf. Eqs. H59I50() ') which gives A^_l(2.5) = —0.06 (cf. Fig. [TB|) . The contribution of the meridian orbits 
has essentially vanished. For rotation parallel to the symmetry axis, only the contribution from the equator orbits is 
modified. For N = 100, a = 0.3, and fuv = 0.6 MeV, the flux $|| = 1.9 (cf. Eqs. |)59I48|) ) which gives A^||(1.9) = -0.3 
(cf. Fig. 1161) . Now the equator orbits contribute with the opposite sign, i. e. valleys become ridges and vice versa. 
The N = 100 mountain, which for w = is generated by the constructive interference of the meridian ridge and the 
equator ridge, becomes a saddle due to the destructive interference. For the same reason, the N — 126 minimum also 
becomes a saddle. The topology does not change for = 50 because the fluxes are smaller. 

In principle, the argument given in the preceding paragraph is not applicable to the region of small a, because 
it is based on the assumption of the two isolated families of equator and meridian orbits. Nevertheless, it describes 
the behavior in this region properly The reason will be given in the appendix, where we shall discuss the region of 
small deformation. It will be shown that the main features of Figs. 1181 can be understood in terms of one family of 
tetragonal orbits, if the integral over the different orientations of these orbits is exactly evaluated. 

Without rotation, the minima and maxima of the shell energy lie on the a = axis, because the spherical symmetry 
correspond to maximal degeneracy of the orbits which is reflected by a maximal amplitude of the oscillating terms. 
Figs. 1171 and 1181 show that for finite rotational frequency the minima and maxima are shifted to slightly negative 
values of a and somewhat smaller values of N. We could not find an explanation for this shift within our simple 
version of POT. 

In the deformed nucleus the shell effects depend on the orientation of the rotational axis with respect to the deformed 
potential. We have studied the most important possibilities that the rotational axis is parallel or perpendicular to the 
symmetry axis. As discussed above, the shell contribution to the rotational energy depends on the flux through the 
orbit. The moments of inertia for rotation about the symmetry axis are determined by the equator orbits and the ones 
for rotation about the an axis perpendicular to the symmetry axis by the meridian orbits. We indicated this in the 
expressions (|52I53II , which we refer to as the parallel and perpendicular moment of inertia. The nucleus rotates about 
the axis with the larger moment of inertia because the rotational energy is smaller. Rotation about the symmetry axis 
corresponds to the generation of angular momentum by sequential particle-hole excitations. The spacings between the 
levels are irregular, and many isomers appear in the yrast line. If the rotational axis is perpendicular to the symmetry 
axis, the yrast levels organize into regularly spaced rotational bands. Hence, rotation selects the contributions of one 
of two types of orbits, which are combined in the ground state shell energy. 

The correlation between the ground state shell energy and the shell contribution to the moment of inertia cause a 
competition between the rotational and the ground-state shell energies, which tend to have the opposite sign. Let us 
discuss two examples of this competition. 

The equator orbits cause the bump (positive shell energy) around N = 106 and a = 0.2 seen in the meridian valley 
(negative shell energy) of Fig. E| This would correspond to a positive Jsh\\ and a negative Jsh±- In this case, the 
difference is large enough to make J'\\ larger than J'± so that the rotation about the symmetry axis is preferred. In 
fact, the N = 106 region is known to be rich in if-isomers, which are states with the angular momentum parallel 
to the symmetry axis. Experimentally, the parallel rotation is not quite as favored as in the calculations. The two 
moments of inertia are about the same, and the yrast line contains high-X bands, for which the rotational axis has 
an intermediate orientation. 

The other example is the region above N — 126. Fig. 1181 demonstrates that oblate shape and parallel rotation is 
energetically favorable as compared to prolate shape and perpendicular rotation (see the right panels, which exagger- 
ate). The reason is as follows. Fig. 1 151 shows that positive contributions to J^sh come from the parallel rotation on the 
oblate side and perpendicular rotation on the prolate side. The J'sh\\ is larger because the factor a|| = 1/2 is larger 
than a± = 1/4 in expression H44|l or, equivalently, the ^-dependence of is stronger than for A4± (cf. Fig. I16|l . 
The macroscopic moments of inertia reinforce this choice. In fact, this region is well known to for high spin isomers, 
which have an oblate shape and the angular momentum aligned with the symmetry axis. Another such region lies 
above = 82. 
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I. Rotation of superdeformed nuclei 

The orbits that cause sheU structure of (prolate) superdeformed nuclei (cf. Sec. lIVPp do not carry rotational flux 
if the rotational axis is perpendicular to the symmetry axis. This is evident for the orbits in the equator plane. For 
the butterfly orbit in Fig. ^Jone must take into account that the rotational flux has a sign. If the rotational axis 
points out of the page, it is positive if the particle runs counterclockwise around the enclosed area and it is negative 
if it runs clockwise. Therefore, the flux has opposite sign for the two wings of the butterfly and the total flux is zero. 
The analogous compensation takes place for the three-dimensional orbits. The projections of the 3D five-point star 
and of the 3D double-traversed triangle on a meridian plane consist of two equal areas contributing with opposite 
sign (cf. Fig. 7 of [sH) cancelling each other. In the observed superdeformed rotational bands the axis of rotation 
is perpendicular to the symmetry axis. Therefore, their moments of inertia should be equal to the rigid body value, 
which is the result of numerous mean-field calculations (see e.g. jT3,E|)- Fig- Elshows that the experimental deviation 
from the rigid-body value is about 5% or less in the A = 150 region. A similar analysis of the A = 190 region is 
not possible, because the moments of inertia are not constant, which likely means the angular momentum where the 
pairing is small has not yet been reached. Other regions of "superdeformed" nuclei have not been included in that 
analysis because the deformation is smaller and it is not clear that the equatorial orbits dominate. Exceptions are for 
the nuclei ^^Tc and ^°®Cd which have an axis ratio close to 2 and indeed have a rigid-body moment of inertia. 

If the axis of rotation is parallel to the symmetry axis the equator orbits carry rotational flux. Since they cause 
the strongly negative shell energy of the superdeformed nuclei, their shell contribution to the moment of inertia will 
also be strongly negative. Therefore, small moments of inertia for rotation about the symmetry axis are expected. 
This means that the appearance of high-_?C isomers and bands is strongly disfavored in superdeformed nuclei. So far 
no evidence for this type of rotation has been reported. 

Concerning the zero-spin shell structure, spherical and superdeformed nuclei are similar. In both cases the periodic 
orbits generate a strongly negative contribution to the shell energy, which generates the shell gap. However, their 
rotational response is different. In spherical nuclei no direction is preferred. There is only one moment of inertia , 
which is roughly proportional to the level density near the Fermi surface and, hence, small. The symmetry axis and 
the short axes play a different role in superdeformed nuclei. The angular momentum has always the direction of the 
short axes, which have the largest moment of inertia . For this orientation, the orbits that cause the shell energy do 
not carry rotational flux. Therefore, the moment of inertia takes the rigid-body value and the shell structure is not 
damped by rotation as in the case of the spherical closed shells. This example demonstrates that the simple relation 
between the level density and the moment of inertia holds only for not too large deformation. In superdeformed 
nuclei, the moment of inertia takes the rigid-body value, although there is a substantial shell gap. 



J. Finite temperature and unresolved spectra 



The study of the unresolved 7-continua permits us to explore regions of finite temperature and large elongation. 
Since the information on the rotational response is less direct than in the resolved spectra, one has to make assumptions 
about the moments of inertia in the evaluation of the data. Therefore, one should know to what extent the moment 
of inertia is expected to deviate from the rigid-body value. 

Finite temperature causes additional damping of the shell structure, which is taken into account by multiplying the 
zero-temperature damping factor D with the temperature damping factor 

Dt = ""^^^/^ (60) 
sinh(7rr^T/?i) ' 

where T is the temperature [H ^ Its argument is nrT/h = TTmLT/{h^kF) = 5.4 x nbTN'^/^ /{2eF) = 

0.38 [MeV]~^A^^^T for the orbit length L = BAR, which determines the basic shell oscillations. The 7-cascades of 
high-spin experiments on nuclei in the mass region 170 are characterized by temperatures of 0.5 MeV or less, which 
corresponds to 0.84 < Dt < 1, i. e. the moments of inertia are not very different from their zero-temperature values. 
The thermal damping becomes important for T = 1 MeV, where Dt = 0.52, and the shell structure is wiped out at 
T = 3 MeV, where Dt = 0.02. These estimates agree well with the microscopic calculations of the shell contribution 
to the moments of inertia at finite temperature in Ref. . 

In the experiments "46] rotational frequencies between 0.4 and 0.8 MeV are reached. For moderate deformation 
(ji = 1.3), the meridian orbits carry rotational fluxes of 1.35 and 2.7, which give modulation factors of 0.6 and -0.1, 
respectively (cf. Figll6|l. Hence, within the frequency interval the shell structure is expected to go to zero and return 
weakly inverted. With increasing angular momentum the nuclei reach the transition to very elongated shapes. Then 
the class of orbits discussed for superdeformation takes over and the moments of inertia take the rigid-body value. 
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The experiments on the unresolved spectra are consistent with weak deviations of the moments of inertia from the 
rigid-body value [46l| . 

K. Currents in the rotating frame 

The deviations of the moments of inertia from the rigid-body value indicate that there must be net currents in the 
body-fixed frame. By the correspondence principle, one expects that these currents are generated by the nucleons on 
the classical periodic orbits. We have not studied the currents in any detail, leaving this interesting question for later. 
In this paper, we give only a heuristic argument. 

For simplicity, let us consider the equator orbits and rotation parallel to the symmetry axis. Classically, the particles 
revolve the symmetry axis, where we consider only square orbits for simplicity (to be brief we omit the subscript □). 
There are two orbits: the particle runs counterclockwise or clockwise, which correspond to Z > and I < 0, respectively. 
The shell contribution to level density associated with the square orbits is g, therefore each of the two contributes g/2. 
Without rotation, there is the same number of particles on both types of orbits, and the net current is zero. When 
we rotate the frame, a particle with angular momentum / gains the energy —luj. Particles will redistribute from the 
unfavored {I < 0) to the favored (/ > 0) orbits, and the net angular momentum and current will no longer be zero. 
For a given chemical potential A, the energy takes a minimum if all particles within the energy interval Ito above A 
flip from —I to I. Hence, the number of flipped particles is lujg/2. Since each gains 21 of angular momentum, the 
total angular momentum gain is l^tog. In agreement with Eq. 155|) . this corresponds to the shell moment of inertia of 

Jsh = gl^. 

This shows that the mass current, and therefore the shell contribution to the moment of inertia, comes from 
nucleons circulating on the contributing classical periodic orbits with the Fermi velocity. Let us consider the current 
distribution. Since the classical mass current of a particle on the square orbit is m/r — pp/ L, the total mass current 
due to all particles on the square orbit is 

Jsh = —lujgsh = ^lujgsh- (61) 
r L 

This current is distributed between the outer circle and the inner circle tangent to the sides of the orbits in Fig. 
[T51 These circles limit the region accessible for a classical orbit. An exponentially decreasing tail will reach into the 
classically forbidden region. 

The current density is j^h = Jsh/'^i where a is the cross section of the area in which jgfi is strong. Perpendicular to 
the symmetry axis the classically accessible area extends between the surface (r = Rg) and r = Rs/y/2 « O.TRg, i.e., 
about 0.3 Rs- For an estimate of extension parallel to the symmetry axis a more profound analysis in the framework 
of POT is needed. As an order of magnitude estimate we take 0.3 i?/. Using the expression 1)21(1 for Ao for the 
total mass current and considering maximal shell contribution (sine function equal to minus one) the current density 
jsh ^ 3.74:mr^'^A^/^uj for rj — 1.3. The current density for the rigid flow is: 

j^i = mpruj - O.^mRsUj- — ^ 0.175TOr~M^/^w, (62) 

where we used r = 0.8i?<;, which is a rough estimate of the average value of r of the particle on the classical orbit. 
The large ratio jsh/jrig 20N~^^^ indicates strong surface currents, which are needed to contribute to the total 
angular momentum an amount comparable to the contribution from the rigid-body flow, which is distributed over 
the whole nucleus. The same type of currents will also be generated by the other polygon-type orbits, both in the 
equator and meridian planes. The shell currents circulate opposite to the rotation if g^h < 0, i. e. they reduce the 
angular momentum. For gsh > they circulate with the rotation and increase the angular momentum. 

V. CONCLUSION 

We have shown that the nuclear moments of inertia at high spins along the yrast line differ substantially from the 
rigid-body value. The differences cannot be attributed to pair correlations, rather they manifest the shell structure. 
Comparing experimental data with quantum mechanical mean-field calculations assuming zero pairing we find a 
similar dependence on the neutron number, which is strongly correlated with the well known shell energy at zero 
angular momentum. The data and the quantal calculations can be interpreted using the semiclassical Periodic Orbit 
Theory, which relates the quantal shell effects to the characteristics of classical periodic orbits in the same potential. 
A number of features, such as the small moments of inertia near closed shells, moments of inertia that exceed the 
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rigid-body value around neutron number 90, the appearance of isomers near neutron number 106, the correlation 
between the deviations of the moments of inertia from the rigid-body values and the ground-state shell energies, and 
the damping of shell structure with increasing angular momentum are explained from this new perspective. The gross 
shell effects persist along the yrast line up to the highest observed angular momenta. The deviation of the mass 
current from the rigid-body flow pattern is generated by nucleons on classical periodic orbits near the nuclear surface. 
The Periodic Orbit Theory provides a qualitative description of these shell effects in terms of classical mass currents 
in the rotating frame. 
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APPENDIX 

Here we demonstrate that the salient features of the shell structure for moderate deformation derive from a simple 
integral over the family of tetragonal orbits. We use the perturbative approach by Creagh 43). Taking into account 
the change of the length of the orbit in linear order of the deformation parameter a (see |4^T|5l|. the contributions 
of deformation and rotation to the action are given by 

A5||/?i = -^kLaP2{cose)a- ^acos0 (A.l) 

AS'i/?i = -ifcLnP2(cos6')a- $nsin6'cos0, (A.2) 

where the Euler angles ip, 9, 4> describe the orientation of the tetragon. As usual, we denote the length of the square 
in the sphere by and the flux through it by The changes of the action due to deformation and rotation give 
rise to the modulation factor 

M = — f [ e''^'^^'^''^'>/^sm9d9d(t), (A.3) 
47r Jo Jo 

which is numerically evaluated. 

Fig. ^]shows the shell contribution to the level density calculated in this way. The similarity with Figs. El El 
1181 is obvious. At large deformation, one may evaluate the integral HA.3I) using the stationary phase approximation. 
The derivative of the Legendre polynomial P2(cos^) is zero for 9 = 7r/2, which corresponds to the meridian orbits, 
and for 9 = 0, tt, which corresponds to the equator orbits. The structure of the level density is the consequence 
of the interference of these two families, as discussed in the main text. However, the interference pattern is still 
recognizable at moderate deformations, where the stationary phase approximation becomes problematic. This justifies 
our interpretation of the shell structure in terms of the interference of the meridian and equator orbits. 
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Ref. [2^ contains a misprint. The factor tt in the argument is missing. 

We ignore the slight difference between a and the deformation parameter e for our qualitative discussion. 
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FIG. 1: Yrast lines for a rotational nucleus ^^"Er (top) and for a non-rotational nucleus ^^^Dy (bottom) together with their 
fit for the ten highest spins. 



FIG. 2: Yrast lines for a selection of nuclei relative to their fit for the ten highest spins: ^^''Er (diamonds), '^^^Hf (squares), 
'^'^^Hf (triangles) and ^^'*Ra (crosses). 



FIG. 3: Experimental moments of inertia as a function of neutron number for small and normal deformation. The symbols 
shown at right indicate the proton number and are the same for all following figures except figure 5. 



FIG. 4: Experimental deviations from the rigid-body moments of inertia as a function of neutron number for small and normal 

deformation. 



FIG. 5: Experimental deviations from the rigid-body moments of inertia as a function of neutron number for superdeformation. 
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FIG. 6: Experimental shell energies at spin as a function of neutron number. 



FIG. 7: Experimental shell energies at spin 20 as a function of neutron number. 



FIG. 8: Maximum angular momentum observed in each nuclei as a function of neutron number. 



FIG. 9: Same as Fig. 0] but for the subset of nuclei with maximum spin greater than 20h. 



FIG. 10; Calculated deviations of the moment of inertia from the rigid-body value as a function of neutron number for different 
shapes (oblate - left panels, prolate - right panels) and different orientations of the rotational axis with respect to the symmetry 
axis (perpendicular - lower panels and parallel - upper panels). 



FIG. 11: Calculated deviations of the moment of inertia from the rigid-body value as a function of neutron number for optimal 
orientation of the rotational axis. 



FIG. 12; Calculated shell energies at spin 30 as a function of neutron number. 



FIG. 13; Classical orbits in the equator plane (upper panel) and the meridian plane of a normally deformed spheroidal cavity. 



FIG. 14; Classical orbits in the equator (upper panel) and the meridian plane of a super deformed spheroidal cavity. 



FIG. 15; Shell energy at zero rotational frequency for two Woods-Saxon potentials. Upper panel; Cavity-like with small 
diffuseness a = 0.05/m, no spin-orbit potential. Lower panel; Realistic parameters from Tab. Qand {N,Z) = (104,78), with 
spin-orbit potential. It is assumed that for each N corresponds Z = inf(78/104). The lines of constant action for the cavity 
are included. The downsloping lines show the rhombi in the meridian plane. The upsloping lines starting at a = —0.2 show 
the squares in the equator plane and the upsloping lines starting at a = 0.2 show the five-point star in the equator plane. The 
dashed lines correspond to maxima and the full lines to minima of the respective shell energies. The values for the action are 
chosen such that for the equator and meridian orbits the lines go through the minima or maxima of the shell energy at a = 0. 
In the upper panel, the action for the five-point star orbit is chosen such that Lt^kp -f = n{2n+ 1/2), n integer, i.e., such that 
the sin-function in (1181 is equal to one. In the lower panel a constant is added, which is chosen such that the magic number 
A'' = 88 for superdeformed shape falls half-way between two lines at a = 0.6. 



FIG. 16; Modulation factors as functions of the rotational flux $. The dashed lines show the quadratic approximation. 
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FIG. 17: Shell energy at finite rotational frequency for a Cavity-like potential rotating perpendicular to the symmetry axis. 
Upper panel: hiu — 0.3 MeV, Lower panel: hu = 0.6 MeV. The lines of constant action are the same as in Fig. 1151 upper 
panel. 



FIG. 18: Shell energy at finite rotational frequency for a realistic Woods-Saxon potential with spin-orbit coupling. Left panels: 
huj = 0.3 MeV, Right panels: hcu — 0.6 MeV, Upper panels: rotation parallel to the symmetry axis, Lower panels: rotation 
perpendicular to the symmetry axis. The lines of constant action are the same as in Fig. llSl lower panel. 



FIG. 19: POT level density (arbitrary units) generated by a the family tetragonal orbits in a rotating spheriodal cavity. Only 
terms of first order in a are taken into account in orbit length. 



TABLES 



TABLE 1: Parameters of the mean-field potential '22\. 





ro(/m) 


a{fm) 




V„ean{MeV) 




Protons 


1.245 


0.599 


0.341 


52.2 


0.746 


Neutrons 


1.260 


0.614 


0.412 
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